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Abstract 

By means of numerical integration we investigate the coherent and incoherent phases in a gen- 
eralized Kuramoto model of phase-coupled oscillators with distance-dependent delay. Preserving 
the topology of a complete graph, we arrange the nodes on a square lattice while introducing finite 
interaction velocity, which gives rise to non-uniform delay. It is found that such delay facilitates 
incoherence and removes reentrant behavior found in models with uniform delay A coupling-delay 
phase diagram is obtained and compared with previous results for uniform delay. 
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I. MODEL SUMMARY 



A. The Kuramoto model 

The popular Kuramoto model of mutual synchronization of coupled oscillators [1] has, 
since its inception, drastically improved the understanding of this prevalent phenomenon. 
Common examples 

as 

include synchronous chirping of crickets, flashing of Chinese fireflies 
4| , clapping of audiences, bursting of neurons, contraction of heart muscles or operation of 
Josephson junction arrays |s ] (s ] to name a few. This model still remains the most succesful 
one, due to its mathematical tractability, combined with the ability to capture the essence 
of synchrony. 

We build up from the definition of the Kuramoto model which is most suitable for direct 
treatment by numerical methods j^ljl]: 

K - 

m = u l + -Y J Mo ] {t)-ut)) (i) 

where i = 1 . . . N, 9i(t) is the phase of the z-th oscillator at time t and u,i are intrinsic 
oscillator frequencies, sampled from yet unspecified probability distribution p(oS) on compact 
support. Kuramoto solved this model exactly in the case of N — > oo and Ui sampled from 
a Lorentz distribution. Solutions for other distributions have subsequently been obtained. 
A model such defined exhibits a (mean-field-type) phase transition between the disordered 
(incoherent) and ordered (coherent) phases as the coupling constant K is increased. Order 
is monitored by the real parameter r, defined as: 

1 N 

r ( t ) e #M = ±^e^« (2) 

When the stationary state is assumed, r(t) = r (re[0, 1]), with r = 1 and r = in total 
coherence and incoherence, respectively. 



B. Introducing delay 

Some real systems cannot be considered without taking delay into account. The popular 
example of a clapping audience synchronizing to clap in unison is valid only for sufficiently 
small audiences, such as opera halls. When distances are of the order of 300m, or higher, 
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the finite speed of sound makes the delay non-negligible. As a result, e.g. football arena 
audiences cannot clap together or have difficulty in coherent singing. 
We start by introducing delay to §1§ in the most general way: 

K N 

The case of uniform delay, = r, is interpreted as coupling of the state at t to the state at 
t — it. The stability of incoherence in such a model has been studied by Yeung and Strogatz 
in 

To introduce non-uniform delay, we arrange the nodes on a square lattice while preserving 
the topology of a complete graph. The coupling remains uniform, however the delay is made 
distance-dependent through the definition: 



1 (sijY 



(4) 



t is interpreted as the inverse velocity. The distance s^- is defined with the so-called "taxi- 
driver's measure" , i.e. as sum of the differences in horizontal and vertical coordinates and is 
measured in number of nodes. To maintain translational invariance we identify the opposite 
edges and the shortest route is always preferred. Hence, when N = L x L, the average 
distance between any pair of nodes (sij) = -|. This definition normalizes the maximum delay 
to 2r and removes dependence on network size. We find that the unmodified parameter r, 
defined by is useful in monitoring the average order in the sample. 



II. SIMULATIONS AND RESULTS 



We have investigated the behavior of the model described by ([3]) and (j4]) by integrating 
the equations ([3]) using a four-step Adams-Bashforth scheme. For simplicity and reference 
with previous results [8|, we set Ui — | for all i, therefore p(u) = 5{u — |). A run for 
one pair of parameters (r, K) consisted of 10000 integration steps with step size At = 0.01, 
out of which the last 6000 were considered for averaging the order parameter r to obtain 
the temporal average r, rejecting the first 4000 when the system is approaching stability. 
Lattices as large as 32 x 32 were considered. The initial conditions, as well as histories of 
9i, were sampled uniformly from [0,27r). 
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FIG. 1. Phase diagram portions for the uniform delay model (a) and the model described by 
(b). The shaded areas visually approximate the incoherent regimes. The phase border points are 
results of numerical integration of (|3|) . 



We have found phase boundaries between the completely ordered and disordered phase 
(Fig. [TJd). For reference, we have produced a diagram for m = r (Fig. [T^,). In the case of 
p(oj) = 8(lj — Uq), the transitions occuring with changing r are instanteous. It is observed 
that non-uniform delay removes the reentrance of synchrony, as intuitively expected. This 
difference is easily understood when considering low coupling. The reentrance in the case of 
uniform delay and K 1 is due to there being sufficiently little difference between states 
at t — t and t for the effect of delay to be approximated by the rotation of all oscillators 
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with average frequency: 



sin(0j(t - r) - 6i{t)) « sin(0j(t) - 4 (t) - cor) (5) 

The reentrance for low coupling then occurs when lot is close to an integer multiple of 2ir 
(r w fcl?, A; = 1,2,...), where the low-delay limit is reproduced. Our simulations reflect this 
heuristic quantitatively up to k = 2 and qualitatively from k = 3 on. In the case of distance- 
dependent delay (Fig. [TJd) and K < 1 the effects of delay are individually approximated 
by phase shifts of ut-^-^ which vary across connections and the low-delay limit cannot be 

\ s ij) 

reproduced by a specific value of delay. Hence, no reentrance occurs. 



III. CONCLUSIONS 



Our results prove that in order to realistically reproduce behavior of synchronizable sys- 
tems in which delay cannot be neglected, the dependence of delay on distance must also 
be accounted for. However, only this dependence was considered in this research, namely 
the case that even distant nodes get to affect other distant nodes with the same strength as 
they affect their closest neighbors, in consistence with the mean-field approximation. This 
is indeed true for some systems, such as digital communication networks. The effect on the 
phase diagram when coupling decreasing with distance is taken into account needs to be 
considered in due course. 
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